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ABSTRACT 

We report a study of the H30a line emission at 1.3 mm from the region around Sgr A* made 
with the Submillimeter Array at a resolution of 2" over a field of 60" (2 parsec) and a velocity 
range of -360 to +345 km s _1 . This field encompasses most of the Galactic center's "minispiral." 
With an isothermal homogeneous HII model, we determined the physical conditions of the ionized 
gas at specific locations in the Northern and Eastern Arms from the H30a line data along with 
Very Large Array data from the H92a line at 3.6 cm and from the radio continuum emission 
at 1.3 cm. The typical electron density and kinetic temperature in the minispiral arms are 3- 
21xl0 4 cm -3 and 5,000-13,000 K, respectively. The H30a and H92q line profiles are broadened 
due to the large velocity shear within and along the beam produced by dynamical motions in the 
strong gravitational field near Sgr A*. We constructed a 3D model of the minispiral using the 
orbital parameters derived under the assumptions that the gas flows are in Keplerian motion. 
The gas in the Eastern Arm appears to collide with the Northern Arm flow in the "Bar" region, 
which is located 0.1-0.2 parsec south of and behind Sgr A*. Finally, a total Lyman continuum 
flux of 3 x 10 50 photons s _1 is inferred from the assumption that the gas is photoionized and the 
ionizing photons for the high-density gas in the minispiral arms are from external sources, which 
is equivalent to ~ 250 09-type zero-age-main-sequence stars. 



Subject headings: Galaxy: center — ISM: individual (Sgr A) 
recombination line — radio lines: ISM 



ISM: kinematics and dynamics — ISM: 



1. INTRODUCTION 

The Galactic center harbors a supermassive 
black ho le (SMBH) with a mass of about 4.2 x 



10 6 M Q (|Ghez et al. ll2008HGillessen et al. Il2009h 



at the position of the radio source Sgr A*. The 
inner parsec of the Galactic center region contains 
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a rich cluster of stars, including at least 55 OB 



stars within the central 0.5 pc (|Paumard et al. 



2006). There are several theories to explain how 
these stars got there. In situ formation, as in gi- 
ant molecular clouds, is improbable because of the 
stars' large velocity dispersion. The accretion se- 
quence may or may not include the formation of an 
intermediate mass black hole in a super-dense stel- 
lar cluster t hat stabilizes the cluster again st tidal 
disruption ([Hansen fc Milosavlievic 1 120031) . The 
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mystery of the formation of these stars is closely 
linked to the balance of mass inflow and outflow 
toward the Galactic center and its interaction with 
the source Sgr A*. Infall of interstellar material 
may lead to star formation up to a certain ra- 
dius from the SMBH, and the subsequent feed- 
back from stellar mass loss may suppress further 
star formation. Both phenomena may determine 
whether Sgr A* is in a quiescent s tate of the active 



galactic nuclei cycle (jLoeb 1 120041 ) 



3) What are the column densities (or emission 
measures), densities, kinetic temperatures, and 
dynamics of the HII gas in the vicinity of Sgr A* 
itself? 

In this paper, we present new observational re- 
sults on the ionized gas in the central parsec of 
the Galaxy based on observations of the H30a 
line made with the Submillimeter Array (SMA0 
and we interpret them in combination with pre- 
vious H92a and continuum measurements made 

Earlier radio recombination line (RRL) studies with the VLA. Section 2 discusses the observa- 

dSchwarz. Breeman. fc van Gorkom II 19891: | Roberts fe Gos&] ons and data reduction. Section 3 presents our 

results and the modeling of the physical conditions 
and dynamics of the ionized gas in the vicinity of 
Sgr A*. Section 4 discusses the ionizing properties 
and possible links to the massive star formation in 
the central parsec, and Section 5 summarizes the 
results and conclusions. We assume a distance of 

At this 

distance, 1 pc corresponds to 24". 



19931 iRoberts. Yusef-Zadeh. fc Goss I Il996h 



re- 
in- 



vealed complex structures of ionized gas 
eluding the "minispiral" and "Bar," whose dy- 
namics remain unclear. Parts of these struc- 
tures have been modeled as gas orbiting Sgr A* 
dSchwarz. Breeman. fc van Gorkom lll989l:ISanders 

1998; Pa umard. Maillard. fc Morris |l2004l:IVollmer fc Dus8hkp c to Sgr A* throughout the paper 
200C ; iLiszt 1 120031 ; iMuzic et al. Il2007tlzhao et al. ' 



2009) while other parts of the structures may be 



gravi tationally unbound ( Yusef-Zadeh. Roberts, fc Birettaj . 
Il998l ). Another possibility is that some of the 
nearby stellar clusters (IRS 13 and IRS 16) have 
powerful winds that, with the help of gravita- 
tiona l focusing, create large-scale, u nbound flows 



OBSERVATIONS 
DUCTION 



AND DATA RE- 



2.1. SMA HSOa Line Data 



(e.g.. lLutz. Krabbe. fc Genzel 111993!) . There are a 



number of limitations in the existing spectral line 
observations of these structures. In the cm-radio 
images made with the Very Large Array (VLA)H 
the recombination line strength is low, the line- 
to-continuum ratio is low, and the bandwidth is 
limited, making high- velocity line wings especially 
difficult to detect. In the infrared, the picture 
given by the emission lines from the ionized gas 
is subject to corrections for dust extinction of ^3 
mag in the central parsec jScoville et al. _ 2003; 



Buchholz. Schodel fc Eckart Il2009t ISchodel et al. 
2010h . A few key questions about the ionized 



medium associated with Sgr A* are: 

1) Is the ionized material related to stars or star 
formation? 

2) Is the ionized material streaming freely in the 
central gravitational potential, or do other mech- 
nisms such as stellar winds and shocks influence 
the dynamics? 



The H30a line observations reported here were 
obtained with the SMA as part of a key project 
on the Galactic center in the period of 2006-2008. 
Our data comprise eleven tracks (Table 1). In 
all cases, the pointing center was the position of 
Sgr A*, and the half-power beam width at the 
231.901 GHz frequency of the H30a recombination 
line was 53" (2 parsec at the distance of Galactic 
center). The 2 GHz bandwidth of the SMA cor- 
relator corresponds to a velocity range of ±1300 
km s . Six tracks were taken in the compact 
array configuration, where the synthesized beam 
was 5.1" x 3.2". Five higher-resolution tracks were 
also acquired, four in the extended configuration 
(1.5" x 1.1" resolution) and one in the very ex- 
tended configuration (0.7" x 0.4"). 

Calibrations of the data were carried out in 



Miriad (jSault. Teuben. fc Wright 111995ft with the 
specific implementation for SMA data reduction!! 



3 The Very Large Array (VLA) is operated by the National 
Radio Astronomy Observatory (NRAO). The NRAO is a 
facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



4 The Submillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is 
funded by the Smithsonian Institution and the Academia 
Sinica. 

5 http: / /www. cfa.harvard.edu/sma/miriad 
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The system temperature correction was made in 
post-processing to compensate for amplitude at- 
tenuation owing to the Earth's atmosphere. Using 
the calibrators listed in Table 1, we determined 
the antenna-based bandpass solutions for each ob- 
serving track, and the bandpass corrections were 
applied to the visibility data. The flux density 
scale was determined from observations of plan- 
ets. Phase corrections were made using the point 
source model of Sgr A*, conveniently located at 
the center of the field, based on self-calibrations. 
The continuum emission was subtracted from the 
spectral data using the Miriad task UVLIN for 
both the point source (~3 Jy) and the extended 
emission (~3 Jy). 

We constructed the H30a line image cube with 
robust weighting (R = 2) (jBriggs Ill995l ). The syn- 
thesized FWHM beam was 1.9" x 1.5" (PA = 26°). 
In order to compare the SMA images with the 
VLA images of the H92a line and 22 GHz contin- 
uum emission, the SMA image cube was convolved 
to a circular beam of 2". The intrinsic resolution 
of the correlator was 3.2 MHz, or 4.2 km s" 1 . We 
smoothed the spectra to a resolution of 15 km s _1 
with a Miriad task. For this paper, we selected the 
velocity range of -360 to +345 km s _1 for analy- 
sis. We found no significant emission beyond this 
range. The rms noise level of each channel image 
was 7 mJy. The image of the integrated H30a 
line is shown in Figure 1 along with spectra in 
the directions of known infrared sources. Figure 
2 shows the distribution of radial velocities deter- 
mined from fitting the peak velocities of the H30a 
line profile in the central 40" of the field. No- 
tice the large velocity gradients along the North- 
ern and Eastern Arms. The Western Arc does not 
show well in the SMA H30a data because most of 
the line emission is below the SMA's 4a sensitivity 
limit. It is detected, however, in single-dish H30a 
data taken with the IRAM 30m telescope. 

2.2. VLA H92a Line Data 

The detailed calibrations and im aging of the 



H92a line data have been described bv lZhao et al 
(|2009l ). They used the line image cube convolved 
to a circular beam (6>fwhm = 2"), the same size as 
that of the H30a line cube. Figure 3 shows a com- 
parison of the spectra of the H92a and H30a lines 
toward selected regions containing IR sources. Ta- 
ble 2 lists the line fluxes and their ratios in these 



spectra. 

2.3. Continuum Data at 22 GHz 

The continuum image at 22 GHz was obtained 
with the VLA high-resolutio n (A and B a rrays) 
data set of 2005 described in IZhao et al. I (|2009h 
and the VLA archival data sets observed in C 
and D arrays during 2004 and 2005. The data re - 
duction procedure described in lZhao et al. (2009) 
was followed. Sgr A* and the X-ray transient 
source (J174540.0290031) 3" south of Sgr A* were 
removed from the visibility data prior to imaging. 
The image of continuum emission was constructed 
from the visibility data with a long-baseline cutoff 
(< 100 kA) to achieve a synthesized FWHM beam 
of 1.95" x 1.41". The dirty image was cleaned with 
the Clark-Steer hybrid algorithm and convolved to 
a FWHM beam of 2", for comparison with other 
data at the same resolution. 

2.4. LRAM H30a Line Data 

We also used additional H30a data taken with 
the HERA multibeam receiver on the IRAM 30m 
telescope, with 11" beams, in order to calculate 
the global parameters for the region. Based on 
the intensity of integrated H30a line emission in a 
lower-angular-resolution (~ 4") image constructed 
by combining the SMA data with the IRAM 30m 
data, we separated Sgr A West into two compo- 
nents: A, the bright minispiral features of the 
Northern and Eastern Arms (/ Sn3Q a dV > 5 Jy 
beam -1 km s _1 ) that trace the high-density ion- 
ized gas, and B, the residuals of Sgr A West after 
subtracting the emission from A. These results are 
in Table 2. 

3. RESULTS 

3.1. Physical Conditions of the Lonized Gas in 
the Central Parsec 

The overall radio properties of the observed re- 
gion (Sgr A West; 80" x 45") determined from the 
H30a, H92a, and 22-GHz continuum images are 
summarized in Table 2. 

3.1.1. Distribution of H30a Line Emission and 
IR Sources 

Figures 1 and 2 show the H30a line data 
in the central parsec. Figure 3 shows a com- 
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parison between the SMA H30a and the VLA 
H92a line spectra in the direction of selected 
IR sources. For optically thin gas in local ther- 
modynamic equilibrium (LTE) with no pressure 
broadening, the ratio of line peak intensity is 



28. 



detailed comparison of the H30a line spectra from 
the SMA and the H92a line spectra observed with 
the VLA yields line ratios in most of the regions 
that are close to the expected LTE ratio, sug- 
gesting that the ionized gas in the central parsec 
is optically thin and under LTE conditions. In 
some regions, however (e.g., IRS 1W, 10W, 16 
and 33), the observed ratios show significant de- 
parture from the expected LTE ratio. 

3.1.2. Distribution of T e 

The equivalent electron temperature, T*, de- 
rived from the ratio of line-to-continuum emis- 
sion under the assumption that the optically thin 
gas is in LTE, can be estimated from our mea- 
surements. In order to avoid contamination from 
both the dust emission, which becomes dominant 
at higher frequencies, and significant synchrotron 
emission at lower frequencies, we used the con- 
tinuum data at 22 GHz, where the contributions 
from both the synchrotron and thermal dust emis- 
sion are negligible relative to the free-free emis- 
sion. Figures 4a and 4b show the integrated 
H30a line image and radio continuum image at 
22 GHz, respectively. We assume that the contin- 
uum emission at 22 GHz is entirely due to opti- 
cally thin free-free emission and scale its intensity 
to 231.9 GHz, according to the power law v~ , 
i .e., a factor of 0.79. Using Equation (1) below 
(jWilson. Rohlfs. fe Huttemeister 1120091) . we deter- 



mined the value of T* , 



J e 

K 



6985 
a(v,T e ) 

Sc 



V GHz/ 1 

1 



\SlAVfwhm/ 



1 



0.87 



N(He)/N(H) 
(1) 



where v is the frequency; Sl is the line flux den- 
sity; Sc is the continuum flux density; AVfwhm 
is the line width (FWHM), and only Doppler 
broadening is involved; N(He)/N(H) is the num- 
ber density ratio, which we take to be 0.08; and 



a(v,T e ) ~ 0.97 (at v = 22 GHz and T e , the elec- 
tron tejn^CTatvrrej^^J^^Klisthe correction to 
the iMezger fc Henderson (1967) power-law ap- 
proximation. The angular distribution of T* is 
shown in Figure 4c. T* varies from 7,000 K to 
10,000 K in most regions of Sgr A West. How- 
ever, in the Bar region (e.g., IRS 2L, IRS 12N, 
IRS 13E, IRS 33), 3" SE of Sgr A*, T* increases 
to ~ 15, 000 K. The significant increase of T* in 
the Bar has al so been derived from li ne observa- 
tions of H92q dRoberts fc Goss Ill993l) an d H66a 
(Schwarz, Brcgm an. fc van Gorkomlll989h . 



3.1.3. Isothermal Homogeneous HII Model 

Since the observed quantities of line and contin- 
uum intensities depend on the density, tempera- 
ture, and path length of HII gas, we can construct 
a model and fit the data to it over a wide frequency 
range in order to determine the parameters of the 
ionized gas. In particular, with the high-resolution 
data obtained from interferometer observations, 
uncertainties due to variations in structure and 
the physical parameters across the source are sub- 
stantially reduced. Thus, the results from such an 
analysis become more reliable. Models for RRL 



ies have been discussed by Puxlev et al. 


(1991); 


Anantharamaiah et al. ( 1993. 200 


J): Zhao et al. 


(1996 


, 1997); Phookun, Anantharamaiah, 


& Goss 


(1998 


) ; and iRodriguez-Rico et al. 


(2004 


)■ They 



found that the main constraints for the models are 
the integrated RRL strength at multiple frequen- 
cies and the observed radio continuum spectrum. 
The observed quantities (line and continuum in- 
tensities) depend in a nonlinear way on the dis- 
tribution of T e and n e , the electron density, both 
along the line of sight and across the telescope 
beam. In the cases of the extragalactic nuclear 
regions, at a given frequency, the line emission of- 
ten arises from HII components under the physical 
conditions that are particularly favorable at this 
particular frequency. Usually, multiple collections 
of HII regions with different physical parameters 
are required to fit the observations. 

Both the SMA and VLA data achieve angu- 
lar resolution (2", ~ 0.1 pc) adequate to re- 
solve the ~ 1 pc minispiral structure. The ef- 
fects due to changes of physical parameters across 
the source are therefore less significant. However, 
the electron temperature and density could possi- 
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bly change along the line of sight; T e and n e de- 
rived from the simple isothermal and homogeneous 
model correspond, therefore, to average physical 
properties of the HII gas in the telescope beam. 
If the HII gas in the beam along the path length 
is isothermal and homogeneous, then the line and 
continuum flux densities from a r egion subtendin g 
a solid angle O are given by (e.g.. LShaver l (|l975h ). 



Sr 



2kv 2 



-VtT P 



T~L I fin + T C 
TL + T C 



and 



(l-e- {TL+Tc ^ -{l-e- TC ) ,(2) 



Ob,, 2 

S c = ^nT e (l~e- TC ) , (3) 



where k is Boltzmann's constant, c is the speed of 
light, and the line (tl) and continuum (rc) optical 
depths are given by 



TL 



Lfv 



pc 



-5/2 



AVh 



km s 



AVpY 
1 + 1.48— M 
AVdJ 



(4) 



and 



TC 



0.08235 



/ n e \ 2 ( Lfv\ t v y 
V cm-V \ pc j ' GHz/ 



cm 

-1.35 

a(i/,T ( 



(5) 



respectively. b n and /3 n are the population depar- 
ture coefficients; n e is the electron density; AVd 
and AVp are the FWHM Dopplcr and pressure 
line widths in km s , respectively; L is the path 
length; and fv is the volume filling factor. Note 
that the last factor in Equation (4) provides a con- 
venient approximation for the transition between 
the pure Doppler and pure pressure-broadened 
cases. 

The Doppler broadening has a thermal and tur- 
bulent component and is given by 



AVn 



'81n(2)fcT e 
M H 



(AV*) 2 



(6) 



where Mh is the mass of hydrogen, and A Vt is the 
turbulent width, which here includes microscopic 
turbulence and dynamical velocity gradients. In 
the central parsec, the thermal width (AT4h = 

J Slnt ^j^'^ -) is less than the turbulent width (AV t ). 
The line width due to pressure broaden ing, com- 
puted bv lBrocklehurst fc Leemanl ( 197ll) . is given 

by 

AVp = 1.7 x 10~ 18 -^7V 7 ' 4 km s _1 , (7) 

J e 

where N is the principal quantum number of the 
transition. The total line width is approximately 
the quadratic sum of AVd and AVp, 



AVt 



FWHM 



V(AV D ) 2 + (AVp) 5 



(8) 



AVfwhm is an observable quantity. We solved 
for the parameters T e , n e , AVt, and Lfv for each 
point in the image in the following way. Since 
AVp -c AVth < AVt for reasonable values of 
T e and n e , we assume initial values of T e and 
n e , and solve for AV t via Equations (6)-(8). We 
computed the values for b n and (3 n for the as- 
sumed values of T K and n e using the NLTE code 

which followed 



of iGordon fc Sorochenkol diOQ 



the analysis of lWalmsley (1990). We then formed 
the flux density ratios i?22GHz = <S , H 30a/5'22GHz, 



^22GHz = S'h92q/'S'22GHz, and R^ia° = Sh30o/Sh92o 

The use of these ratios eliminates the solid angle 
CI in Equations (2) and (3). At each point in the 
image, we used these ratios to solve for T e , n e , and 
Lfv We then iterated the process, i.e., solved for 
a new estimate of AV t with these values of T e and 
n e . The process converged in a few iterations. 

The parameters derived for the line of sight 
for each of the regions are summarized in Table 
3. Instead of Lfv, we report the emission mea- 
sure {EM = n\Lf v ). The ratio T e /T* is given 
in column 6, which mainly corrects for the ef- 
fect due to underpopulated electrons at the level 
N = 30 as compared to that expected under the 
LTE; the ratio AVt /AVfwhm in column 7; and 
AV P H92Q /AV D rH , the ratio of the H92a line width 
of pressure broadening to that of thermal Doppler 
broadening, in column 8. The pressure broadening 
term (AVj? 30 ") is negligible. 

The uncertainties in the physical parameters 
mainly come from the uncertainties in the mea- 
surements. Although the relation between the 
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physical parameters and observable quantities is 
nonlinear, we assess the errors through linear ap- 
proximation. We use the optically thin assump- 
tion, which appears to be valid at the frequencies 
used in this paper. The fractional uncertainty in 
T e , which is proportional to the fractional uncer- 
tainties of the measured ratios, i?22GHz an d ^h92qj 
lies in the range from 7% for IRS 1W to 39% for 
IRS 12N. The fractional uncertainty in n e is dom- 
inated by those of observed line ratios -Rh92q an d 
ranges from 6% for IRS 1W to 36% for IRS 12N. 

The variations of the H30a line intensities re- 
flect the changes of excitation conditions in the 
region. Our model fitting suggests that, given sim- 
ilar solutions in n e and T e , the large difference in 
the H30a line intensity between the two extreme 
regions (IRS 1W and IRS 7) is mainly due to the 
difference of the emission measure (EM) or the 
path length along the line of sight (see Table 3). 

In short, the line ratio of H30a to H92a varies 
from 15±4 to 32±4. The derived mean values of 
the line optical depth are Ti(H30a) = —0.001 
and tl (H92a) = —0.006 with standard devia- 
tions of ct Tl (H30q) = 0.0004 and cf TL (W2a) = 0.005 
from the corresponding emission regions. The 
free-free optical depth at the H30a line frequency 
t c (231.9GHz) varies from lxl0~ 4 to lxl0~ 5 
while t c {8.3 GHz) is in the range between 0.16 
and 0.01. Because of the small optical depths in 
both line and continuum, the non LTE effects in 
the minispiral, in general, are not critical. 

There is a diffuse nonthermal sync hrotron com- 



pone nt in the central parsec region (Ekcrs ct al 



19831 ). Based on the flux dens ity of 22 Jy a t 20 cm 



in the central 40" x 70" from lEkers et al. I ([1983t h 



we estimated that the contribution to the contin- 
uum flux density in a 2" beam is ~10 mJy under 
the assumption that the nonthermal spectral index 
a ~ —0.5 (Sv oc v a ) and the emission are evenly 
distributed. The optical depth of the H92a line 
is small (tz, ~ —0.01). The contribution from the 
stimulated H92a emission of ~ 0.1 mJy is insignif- 
icant in the regions discussed in this paper. The 
contribution from the stimulated line emission by 
the dust emission at 232 GHz is negligible because 
of the small optical depth of the H30a line. 

The observed line widths of the HII gas in the 
minispiral arms at the Galactic center are dom- 
inated by AVt, which is consistent with velocity 
gradients due to gas motion around Sgr A* or the 



SMBH being primarily responsible for the broad- 
ening in the line profiles observed in the 2" beam. 

3.1.4- Northern and Eastern Arms 

The ionized gas of the high-density compo- 
nents in the Northern and Eastern Arms has a 
mean kinetic temperature of T e k 6 x 10 3 K 
and an electron density of n e in the range of 



10 



4-5 



cm , which is consistent with the values 



(10 4 cm 3 ) derived from Paa line emission at 1.87 



lim (jScoville et al. 2003) for a given uncertainty 
in extinction corrections in the near IR. The HII 
gas properties in the Northern and Eastern Arms 
are similar to these in HII regions around young 
and/or evolved massive stars. 

3.1.5. The Bar 

We confirm that the ionized gas in the Bar 
(e.g., IRS 2L, IRS 12N, IRS 13E, and IRS 33 in 
Table 3) appears to have a higher kinetic tem- 
perature in the range of 8.5 x 10 3 K < T e < 
1.3 x 10 4 K, as suggested by e arlier RRL observa- 
tions ([Roberts fc Goss Ill993l ). The electron den- 
sity shows a range from a few times 10 4 to a few 
times 10 5 cm~ 3 . 

The heating process may be complicated in the 
Bar because of complex gas interactions. There 
are several energy sources in the region. First, 
strong winds from the central stellar clusters such 
as IRS 16 and IRS 13 could provide kinetic energy 
to compress the ionized gas in this region. The 
heating due to the interaction between the winds 
and the ionized gas in the Bar likely occurs in the 
regions localized to the cluster sources. 

In addition, the well-known minicavity (marked 
with the gray circle in Figures 4a and 4b) , a region 
depressed o f radio continuum, is located 3.5" SW 
of Sgr A* (|Yusef-Zadeh. Morris, fc Ekers I Il989h . 
This region is likely created by a fast wind outflow 
from the sources in the vicinity sweeping up the 
gas in the orbiting ionized streams, as suggested 
by the morphology of a bright [Felll] lin e in the in- 
frared (|Lutz. Krabbe, fc Genzel 1993 ). However, 
the minicavity is a few arcsec displaced from the 
high-kinetic-temperature region, which in fact is 
best seen as a depression associated with H30a line 
emission. The minicavity can also be explained by 
Melia. Coker. fc Yusef-Zadeh (|l996h as a down- 
stream, focuccd flow due to the accretion by the 
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SMBH at Sgr A*. 

Alternatively, the high temperature in the Bar 
could result from a collision between two ion- 
ized orbiting flows, as suggeste d in the dynami- 



cal model of Zhao et al. ( 2009t ). A computation 



using their Keplerian orbital model (see Section 
3.2) shows that the specific kinetic energy reaches 
a maximum (see Figure 4d) in the Bar region 
where high kinetic temperatures (see Figure 4c) 
are observed. The distribution of specific kinetic 
energy (Figure 4d) shows a substantial increase 
from 2 x 10 14 erg g _1 at a region near IRS 5 in 
the Northern Arm to 12 x 10 14 erg g" 1 at IRS 33 
in the Bar. However, we note that a high kinetic 
energy of the orbiting flow does not have to corre- 
spond to a high kinetic temperature; the ionized 
flows can release their kinetic energy to shocks if 
the ionized flows collide with each other. If they 
collide with the strong stellar winds from the clus- 
ters of stars in the Bar area, then the gas could 
be heated by the shocks to higher temperatures. 
In addition, the model also shows that the ionized 
flows of the Northern and Eastern Arms do col- 
lide at the region with high kinetic t emperatures 
(see Figure 21d of IZhao et al~l (|2009h and the 3D 
structure of the minispiral arms in Figure 5). The 
observed high kinetic temperatures in the Bar sug- 
gest that the ionized gas could be further heated 
by the shocks generated from the collision. 

3.2. Kinematics and Dynamics 

Kinematics of the ionized gas in both the 
Northern Arm and Western Arc have been ex- 
plained as Ke plerian circular motion of a circum 
nuclear ring (Schwarz, Brcgm an. fc van Gorkom 



other kinematic features, the "Western Bridge" 
and the "Bar Overlay." The VLA observations of 
the H92a line suggest that the loci of the three 
minispiral arms can be fitted with three families 
of elliptical orbits; the orbits of both the North- 
ern and Eastern Arms have high eccentricities 
while the Western Arc is in nearly circular mo- 
tion (Zh ao et al. 20091 ). The VLA measurements 
of the proper motions and the radial velocity re- 
solved the degeneracies in the inclination angles, 
suggesting that the two ionized flows from the 
Northern and Eastern Arms collide in the Bar 
region. 



On the other hand, iMuzic et al. I (|2007l ) sug- 
gested that the non-Keplerian motions might be 
important in the kinematics of the filaments ob- 



served near Sgr A*. In addition, IVollmer fe Duschl 



(|2000h noticed that a substantial amount of veloc- 
ity dispersion along with additional inwards radial 
motion is required in order to fit the observed ve- 
locities of the ionized flows with their Keplerian 
circular-orbit model. 

3.2.1. 3D View of the Minispiral Structure and 
the Circumnuclear Disk 



Usi ng the Keplerian model described in lZhao et al. 
(2009) along with the orbital parameters derived 
from fitting the VLA H92a data and proper mo- 
tion measurements, we recomputed the 3D loca- 
tions of the three ionized arms. Figure 5 shows 
a 3D visualization of the three ionized flows in 
the Northern, Eastern, and Western Arms (Arc), 
which illustrates that the Northern Arm and West- 
ern Arc are nearly coplanar based on the similar 
long itudes of their ascending nodes (f2) and in- 



1989; lLacv. Achtermann. fe Serabvn 1991; IVollmer fe Duschha tion angles (i) derived from fitting the loci 



200oi lLisztl l2003h. Based on the VL A observa 



tions of H92a. lRoberts fe Goss I (|1993I ) found that 
the ionized ring consists of three separated kine- 
matic components, the Western Arc, Northern 
Arm, and Bar. The radial velocities of the West- 
ern Arc can be fitted with a Keplerian model with 
circular motion. Using the observations of the 



Brack ett-7 line, iPaumard. Maillard. fe Morris 



of the ionized flows. The mean values of the two 
orbital angles, U = 71° ± 6° and i = 119° ± 3°, 
are inferred for the two coplanar orbits (North- 
ern Arm and Western Arc), which form the ion- 
ized ring. The orbital plane of the ionized flows 
appears to be aligned with the pl ane of the cir 
cumnu clear disk (CND) derived by lGiisten et al. 



(|2004l ) found that the kinematics of the ionized 
gas in the Northern Arm can be explained as a 
bundle of elliptical orbits. In addition, the Eastern 
Arm may be interrupted by a high-velocity com- 
ponent "tip". Also, the Bar component extends 
to the IRS6/34 region which is affected by two 



(|1987I ) from molecular line observations. If the 
fitted parameters fr om both the molecular lines 
(IGiisten et al. Ill987l) and hydrogen recombination 
lines (Schwarz. Brcg man. fe van Gorkom I Il989t 



Lacv. Achtermann. & Serabvn 1991; Roberts & Goss 
19931: IVollmer fc Duschl l2000HLiszt ll2003l : IZhao et al. 



2009) are correct, the Keplerian motions dominate 
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the kinematics in the central parsec, and the plane 
of the main nuclear gaseous disk or the CND can 
be described by the mean orbital angles inferred 
above. Therefore, the orbital plane of the East- 
ern Arm (i = 122° ± 5° and Q = -42° ± 11°) is 
nearly perpendicular to that of the CND. Figure 5 
shows that the Eastern Arm flow runs into the 
CND from the northeast and the near side with 
respect to Sgr A* and collides with the Northern 
Arm flow in the region located south and behind 
Sgr A*. 

3.2.2. Kinematics and Keplerian Motions 

Figure 6 shows the spectra taken from 12 po- 
sitions along each of the two ionized arms (shown 
in Figure 7). The observed spectra are com- 
pared with the radial velocities (green bars) de- 
rived from the Keplerian model. Note that we 
assume that the LSR velocity of the black hole, 
V z . is to within ± 5 km s _1 (see discussion 



Gillessen et al. 



(l2009h . In general, the kine- 
matics of the ionized gas follow the Keplerian or- 
bital motion expected for the gravity of the SMBH 
(4.2 x 10 6 M m ) at Sgr A*. Models of the enclosed 
mass (see Oh. Kim, fc Figer 2009t ) show that the 



added enclosed mass starts to rise above the black 
hole's mass starting at a radius of 0.3 pc (7") from 
Sgr A*. At a radius of 1 pc (24") from Sgr A*, 
near the tip of the minispiral seen in projection, 
the enclosed mass has nearly doubled, to about 
8 x 10 6 Mq. Thus, for a given circular orbit with 
a radius of 1 pc the actual orbital velocity would 
be about 1.4 times higher than in a pure Keple- 
rian approximation based on the mass of the black 
hole only. Apart from this effect, however, in sev- 
eral regions the radial velocities deviate even more 
significantly from those predicted from the Keple- 
rian model. 

3.2.3. Non-Keplerian Components 

In the northwest end of the Eastern Arm near 
the location E17 (Figure 7a), a feature with radial 
velocities blue-shifted up to ~ —200 km s _1 shows 
a large deviation from the value of ~ —100 km s _1 
predic ted by the Kep l erian model (Figure 6). 
From IPaumard et al. I (|2006l) . we find that at 
least five massive stars— IRS 34W (Ofpe/WN9), 
IRS 34E (09-9.51), IRS 34NW (WN7), IRS 7SW 
(WN8), and IRS 3E (WC5/6)— are located in a 
region with a radius of 2.2" (indicated with a 



large black circle near region E16 in Figure 7a), 
referred to as IRS 34 hereafter. Except for IRS 
3E, which has no proper motion measurements, 
the p roper motion data of the remaining four 
stars ( Paumard et al. 20061) show that the stars 
all move southwest toward the ionized flow with 
a mean transverse velocity of 185 km s _1 (PA 
— -134°) (as shown by an arrow in Figure 7a). 
We compare the radial velocity of the stars with 
the spectrum (inset c in Figure 7a) of the H30a 
line emission integrated from the IRS 34 region. 
The four stars besides IRS 3E show negative ra- 
dial velocities in the range -150 km s _1 to -340 
km s _1 , which is blue-shifted with respect to the 
peak velocity (—146 ± 4 km s _1 ) of the broad 
(AVpwhm = 160 ± 8 km s _1 ) spectral feature. 
The stellar winds from the Wolf-Rayet (WR) stars, 
evolved from early type O-stars with mass loss 
rate of 10~ 5 M Q yr _1 , appear to play a consid- 
erable role in the local kinematics. The strong 
stellar winds with velocities up to a few times 10 3 
km s _1 could sweep up the ionized flow of the 
Eastern Arm, causing the peculiar velocities and 
the velocity gradient across the Eastern Arm at 
the northwest end. In fact, the motions of the 
filaments associated with the minispiral found by 
( Muzic et al. 1120071 ) in infrared also suggest that 
a model of purely Keplerian motions of the gas is 
inadequate and the filaments could be due to the 
interaction of a fast wind with the minispiral. 



In the tip region ([Paumard. Maillard. fc Morris 

2004) with peculiar radial velocity that deviates 



signficantly from the Keplarian velocity by +150 
km s _1 (see panel El 2 in Figure 6 and Figures 
7a and 7b), there is also a WN8 star, IRS 9W. 
The radial velocity (+140 km s" 1 ) of IRS 9W 



( Paumard. Maillard. fc Morris I 120041 ) is similar 



to that of ionized gas in the Eastern Arm. A com- 
parison between the stellar transverse velocity of 
215 km s" 1 (PA = 51°) and that of 304 km s" 1 
(PA = -67°) predicted from the Keplerian model 
for the ionized flow at the location of IRS 9W 
suggests that the star does not follow the orbit of 
the Eastern Arm but may run across the Eastern 
Arm. Then, the high-radial velocity tip is likely 
to be the part of the gas in the ionized flow that 
is compressed and accelerated by the suspected 
strong wind from the IRS 9W star. The velocity 
discontinuity at the tip may be evidence for the 
presence of shocks in the region. 



S 



In addition, deviations from Keplerian motions 
are also observed in the Bar region. The non- 
Keplerian kinematics could be attributed both 
to the interaction between the ionized flows in the 
Northern and Eastern Arms and to the interaction 
of the ionized flows with the strong winds from the 
massive star clusters (the IRS 16 and IRS 13 clus- 
ters, for example). Well-organized polarized emis- 
sion from magnetically aligned dus t grains in the 



the confluence of gravita tionally focussed winds 



central par s ec has been observed (lAitken et al 



199lL 119981: iGlasse. Aitken fc Roche I l2003h 



sug- 
gesting a magnetic field strength of >2 mG. The 
corresponding energy densities in the magnetic 
fields are ^ 1.6 x 10 -7 erg cm -3 , about twice larger 
than the thermal energy densities in the ionized 
gas (T e ~ 7 x 10 3 K and n e ~ 6 x 10 4 cm" 3 at IRS 
1W) but two orders of magnitude less than the ki- 
netic energy densities of the orbital motions. The 
magnetic fields frozen in the ionized streamers can 
not substantially alter the bulk motions of the gas 
but produce some cumulative magnetohydrody- 
namic (MHD) effects, such as orbital compression 
of the streamers as implied from the observed ev- 
ide nce for the convergence of field lines near Sgr 
A* (|Glasse. Aitken fe Roche"ll2003l ) . and a helix of 
a natural morphology for a twisted magnetic field 
frozen into a plasma as observed in the North- 



ern 



. Arm (lYusef-Zadeh. Morris, fc Ekers I 119891: 
Zhao et al. Il2009h . In the region close to Sgr A*, 
the magnetic field probably becomes stronger, and 
the MHD effects may also play a significant role 
in the kinematics of the ionized flows. 

An interesting anomalous feature located ~ 
2" northwest of Sgr A* can be seen in the inte- 
grated H30a line intensity image (Figure 1). The 
radial- velocity image (Figure 2) indicates a large 
velocity gradient across this feature. The H30a 
line spectrum (inset d in Figure 7) suggests that 
this feature is composed of at least three spec- 
tral components (> 4<r) at radial velocities of 160 
km s _1 , 39 km s^ 1 , and -250 km s _1 , with a line 
width of ~ 50 km s _1 , a typical velocity width of 
the ionized gas in the minispiral arms. We note 
that there is a compact contin uum source called e 
(|Yusef-Zadeh fc Morrislll99ll ). located 1" south of 
this anomalous feature, corresponding to a hole in 
the H30a line images (e.g., Figure 7a and Figure 
1). The source e appears to be moving away from 
Sgr A* with a tra nsverse velocity of 340 km s _1 
(jZhao et al. Il2009h . which has been interpreted as 



from the IRS 16 cluster (Wardle & Yuscf-Zadeh 



Il992h . The H30a line property of the anomalous 
feature appears to be consistent with the focused- 
wind model, suggesting that the velocity compo- 
nents in the anomalous feature correspond to ion- 
ized blobs of the minispiral arms that have been 
blown out of the region e due to the high pres- 
sure of the confluent stellar winds. However, given 
the closeness to Sgr A*, the detection of this fea- 
ture needs to be further confirmed with higher- 
sensitivity observations. 

Finally, in the locations with large angular off- 
sets from Sgr A*, we noticed that the deviations in 
radial velocity from the Keplerian motions become 
significant for the Northern Arm (> 50 km s~ 4 ), 
as shown in both Figures 6 and 7, while for the 
Eastern Arm, the deviations are less significant. 
The difference in the deviation of the observed 
radial velocities from the Keplerian velocities be- 
tween the Northern and Eastern Arms suggests 
that the distribution of mass has a preference for 
the plane of the CND. 

4. IONIZATION 

4.1. The High-Density Gas in the Northern and 
Eastern Arms 

We have shown that the H30a lines observed 
with the SMA are excellent tracers of the high- 
density HII components in the Northern and East- 
ern Arms. Based on the observations and derived 
physical parameters, we can estimate the ionizing 
photon rate (Q = nlasV) required for maintain- 
ing the ionizations in these regions for the case 
where the ionized gas is optically thick in the Ly- 
man lines and where as is the ionization coef- 
ficient, and V w VlD 2 Lfv is the volume of the 
ionized clumps. The ionizing photon rate can be 
written as 



Q{H) = 1.31 x 10 4 



afl(T e ) 



'B 



arcsec 
Lfv 



a B (10 4 K 



10- 3 pc 



10 4 cm" 3 



phot, s 



D 



8 kpc 



(9) 
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where 6b is the resolution of the observations 
(beam) with the assumption that source fills 
the beam; and a_B(T e ), the total recombina- 
tion coefficient to excite d levels, is given by 
Hummer fc Seatoiil ( 19631 ). We calculated the val- 
ues of the inferred ionizing photon flux required to 
maintain the ionization of the high-density com- 
ponents in the Northern and Eastern Arms as 
well as the HII mass (Mhii) and excitation pa- 
rameter (U). On the assumption that the HII gas 
is internally ionized, the equivalent zero-age main 
sequence (ZAMS) star that is the ionizing source 
was estimated. The internal-ionization model re- 
quires at least ~ 25 09-type stars embedded and 
roughly evenly distributed in both the Northern 
and Eastern Arms. 

4.2. Ionizing Sources 

We now examine where the known OB stars lie 
in the region. Figure 8 shows the 90 massive stars 
within a radius of 14" (0.5 pc ) of Sgr A* identi- 
fied by IPaumard et"al~l (|2006h . among which 55 
are OB stars. The 14 OB stars within a radius of 
0.85" (0.03 pc) from Sgr A* are known as S stars. 
No obvious correlations appear between the loca- 
tion of 41 OB stars outside the radius of 0.85" and 
the distribution of the high-density HII gas in the 
minispiral arms on the projected plane. If we as- 
sume the average solid angle for the minispiral 
subtended to be f2 s , then the required number of 
ionizing photons from the stars will be 



Q*(H) = Q(H){4ic/n s 



(10) 



If we assume 4ir/il s ~ 10, then 250 09-type stars 
are required. 

For the lower-density (n e ~ 8x 10 2 cm -3 ) com- 
ponent B, the ionized mass of ~ 340 Mq is in- 
ferred, requiring a Lyman photon flux of 1.2 x 10 50 
phot, s^ 1 or -100 09-type stars (ZAMS). The 
250 09-type stars required for externally main- 
taining the ionization of the minispiral appear to 
be adequate in supplying the leaking Lyman pho- 
tons from the minispiral to further ionize the lower 
density component that pervasively distributes in 
the region. Thus, the ionizing source to maintain 
the ionization state of Sgr A West is equivalent 
to -250 09-type stars (ZAMS). Considering that 
a substantial amount of the ionizing photons can 



be produced by the hot evolved ma ssive stars and 
the survey of Paumard et al. ( 20061 ) is incomplete, 
there may be enough stars to ionize the region. 

The origin of the ionized minispiral arms and 
the birthplace of the massive stars found in 
the central parsec are still matters of conjec- 
ture. There are two possible models that may 
explain the presence of the massive stars near 
the SMBH. First, the stars formed elsewhere 
outside the cen tral parsec and migr ated to the 



Galactic center ( Gerhard et al. 1 120011 ). However 



compared with the age of the massive stars in 
the central parsec, it takes too long for stars 
to move to the center from several parsec away 
by losing their orbital energy through dynamical 



fricti on (Portcgies Z wart. McMillan, fe Gerhard 



20031 : iKim fc Morris 1 120031 ). Alternatively, the 
stars may have formed near the SMBH by an un- 
usual astrophysical process, such as the fragmen- 
tation of a massive disk rotating near t he SMBH 



([Navakshin. Guadra. fc Springel 1 120071). Recent 



modeling and numerical simulation s (Navaks hin. Guadra. fc Springel 



2007HWardle fc Yusef-Zadeh l2008tlBonnell fc Rice 



12008 ) show that for an infalling giant molecu- 
lar cloud (10 4-5 M ) interacting with an SMBH 
(1-3 x 10 6 Mq), a small fraction of the gas cloud 
becomes bounded to the black hole, forming an ec- 
centric disk that quickly fragments to form stars. 
The tidal disruption of the self-graviting clouds 
formed coherent spiral structures, i.e., the mini- 
spiral arms, which are ionized by the Lyman con- 
tinuum photons produced from the newly formed 
massive stars. 

5. SUMMARY AND CONCLUSION 

We report results from the SMA observations of 
the H30a line emission from the Galactic center at 
a resolution of 2". We measured the velocity pro- 
files over the central 60" (2 parsec) region around 
Sgr A*. We carried out an analysis of the line and 
continuum data based on an isothermal homoge- 
neous HII model for each component. With the 
constraints from the SMA H30a data at 1.3 mm 
and the VLA H92a data at 3.6 cm and continuum 
data at 1.3 cm, we determined the physical condi- 
tions of ionized gas in the Northern and Eastern 
Arms of the minispiral. We found that the typi- 
cal electron density and kinetic temperature in the 
arms are 3-21 x 10 4 cm" 3 and 5,000-13,000 K, re- 



10 



spectively. The highest density of 2.1 x 10 5 cm -3 
occurs in the IRS 13 region. Higher temperatures 
up to 13,000 K are found in the Bar region. Both 
the H30a and H92a line profiles are broadened 
due to large velocity gradients along the line of 
sight and across the 2" beam produced by the 
dynamical motions near Sgr A*. For the H92a 
line, the line width due to pressure broadening ap- 
pears to be comparable to that of thermal Doppler 
broadening, while for the H30a lines, the pressure 
broadening is negligible. 

Using the orbital parameters derived under the 
assumption that the three minispiral flows are in 
Keplerian motion, we calculated the 3D geometry 
of the minispiral structure. We showed that the 
ionized flows of Northern Arm and Western Arc 
are almost coplanar, and the plane of the ionized 
flow in the Eastern Arm is nearly perpendicular 
to that of the Northern Arm and Western Arc. 
The Eastern Arm flow collides with the Northern 
Arm flow in the Bar region located 0.1-0.2 parsec 
south of and behind Sgr A*. The high kinetic 
temperature gas in the Bar is probably due to the 
heat in shocks generated by collisions. 

We compared observed radial velocities in the 
Northern and Eastern Arms with those computed 
from the Keplerian model with a mass of 4.2 x 
10 6 M e of the SMBH at the position of Sgr A*. 
We show that Keplerian motion dominates the 
kinematics of the minispiral arms. In some re- 
gions, significant deviations exist between the ob- 
served radial velocities and those computed from 
the Keplerian model. Near the IRS 9W and IRS 
34 regions, the observed peculiar velocities in the 
minispiral are likely a result of the ionized gas 
being swept up by the strong stellar winds from 
nearby WR stars. 

A total Lyman continuum flux of 3 x 10 50 



phot. 



is needed to maintain the ionization of 



the gas in Sgr A West with an inferred ionized 
mass of 350 M , which requires ~250 09-type 
(ZAMS) stars. About 10% of the total number 
of UV photons is required to externally maintain 
the ionization of the high-density gas in the North- 
ern and Eastern Arms, which have an ionized mass 
of 12 Mc^, less than 5% of the total ionized mass 



(jLiszt Il2003h 
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Fig. 1. — The image of the H30a line intensity integrated from the velocity range between —360 and +345 
km s _1 with the cutoff of 30 mJy beam" 1 (~ 4a) in each of the channels. The color wedge shows the intensity 
in the range between to 10 Jy beam -1 km s _1 . The beam size (#fwhm = 2") is shown in the bottom 
right. The dot-cross at the coordinate origin marks the position of Sgr A*, where («j20oo = 17 h 45 m 40.037 s , 
<5j20oo = — 29°00'28.11") defines the origin of the coordinate system. The plus signs mark the positions from 
which the H30a line spectra in the beam are taken and a few additional sources discussed in the text. The 
red curves are the Gaussian fits (multiple components in some cases) obtained from a least-squares analysis. 
The spectra are labeled with the names of their corresponding IRS sources. 
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RA offset (arcsec; J2000) 

Fig. 2. — The distribution of the radial velocities determined from fitting the peak velocity of the H30a line 
image cube observed with the SMA in the central 40" region. The color scale denotes velocity from -360 to 
+345 km s" 1 . The circular beam #fwhm = 2" is shown at the bottom right. The dot-cross at the coordinate 
origin marks the position of Sgr A*. 
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Fig. 3. — A comparison of the H30a line spectra (red) with the VLA H92a line spectra (green) taken from the 
selected IRS source regions with a resolution of #fwhm = 2". The VLA H92a spectra have been multipled 
by a factor of ^h30c</i / H92q = 28. The vertical bars indicate the uncertainty of 3c. 
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Fig. 4. — a) Integrated H30a line emission observed at 1.3 mm with the SMA. b) Radio continuum image 
made from VLA observations at 1.3 cm with a resolution of #fwhm = 2". The radio emission associated 
with Sgr A* and the X-ray transient (J174540. 0290031) has been removed, but the position of Sgr A* is 
marked by the dot-cross at the coordinate origin in each of the figures, c) Pseudo-color shows the LTE 
temperature (T*) derived from images of H30a line emission and radio continuum based on Equation 1. 
d) The image of the s pecific kinetic energy computed using the Keplerian orbital model described in Table 5 
of IZhao et al~l (l2009h . The color scale shows the specific kinetic energy in units of 10 14 erg g 1 . In Figures 
3c and 3d, the electron temperature (T e in 10 3 K) and the electron density (n e in 10 3 cm' 3 ) are shown 
at the positions of selected IR sources. These numbers are derived using the isothermal homogeneous HII 
model with the non-LTE effects. The results are tab ulated in Table 3. The thick gray cir cles in (a),(b), (c) 
and (d) indicate the well-known minicavity region (jYusef-Zadeh. Morris, fc Ekers 1 119891 1: the black circles 
in (c) and (d) indicate high-temperature regions. 
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Fig. 5. — 3D view of the minispiral structure and velocity along each of the minispiral arms in the nucleus 
of the Galaxy center ed at Sgr A* (black dot), which are computed using the Keplerian model described 



Zhao et al. 



_ ( 20091) . T he orbital pa r amete rs were derived from fitting the VLA H92a line data and are 
summarized in Table 5 of IZhao et al. ] (<2009l) . The arrows indicate the direction of the ionized flows along 
the orbits. The color illustrates the radial velocity between —250 to +250 km s^ 1 (color wedge at top). The 
gray shadow indicates the loci of the ionized flows projected on the sky plane. 
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Fig. 6. — Left: Spectra of the H30a line taken along the Northern Arm. Right: Spectra of the H30a line 
taken along the Eastern Arm. The locations where the spectra were taken are shown in Figure 7. The green 
bars mark the radial velocities along the Nort hern and Eas t ern A rms computed using the Keplerian model 
with the orbital parameters given in Table 5 of lZhao et al. (2009), assuming that the mass of the SMBH is 
4.2 x 10 6 M & and its radial velocity V z = km s _1 with respect to the LSR. The thin blue bars (in panels 
E12, E16, and E17) indicate peculiar velocities possible due to interactions with the strong stellar winds 
from the nearby WR stars (see Section 3.2.3). Note that the spectrum marked "E13" and shown in red in 
both columns comes from a common position, with the appropriate Keplerian velocity marked for each arm. 
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Fig. 7. — A comparison between the observed radial velocities (a) and those computed from the Keplerian 
model (b). The black circles along the two spiral arms mark the locations of the H30a spectra displayed in 
Figure 6. The red circle indicates the high temperature region. The dot-cross marks the position of Sgr A* 
at the coordinate origin. The inset (c) is the line spectrum of H30a toward the IRS 34 region indicated with 
a large black circle, which includes the five massive sta rs IRS 7SW (WN8), IR S 34W (Ofpe/WN9), IRS 34E 
(09-9.51), IRS 34NW (WN7), and IRS 3E (WC5/6) dPaumard et al. Il200fih . The mean proper motion of 
the first four stars (and a fifth star without proper motion) is indicated with an arrow near the E16 region 
in (a). The radial velocit ies of the five st a rs are marked with bars in inset (c). The proper motion of the 
IRS 9W star (WN8) from lPaumard et al. 1 (120061) is indicated with an arrow near the E12 region in (a). The 
inset (d) shows the H30a line spectrum toward an anomalous feature (—1.5", 0.5") northwest of Sgr A* and 
1" north of e. 19 



o 

C<1 




o 
o 

CD 
0.2 

!-5 



O 



o 

CD 
m 
O 

U 



o 



CL) 
t» 

*H — I 

o 
o 

H 
Q 



O 



o 

C\2 



20 



10 -10 

RA offset (arcsec; J200Q) 



-20 







5 



10 



15 



20 



25 



Fig. 8. — The 90 massive stars of Paumard et al. ( 20061 ) are overlaid on the H30a line image. Among the 
90 massive stars, 14 OB stars (orange), known as S stars, are located within a radius of 0.85" from Sgr 
A* (white plus sign), 41 OB stars (white) are outside the S cluster, and the rest are Ofpe/WN9 WR stars 
(green). The cross- and dot-circles represent the stars with proper motion data, and the open circles indicate 
stars with no proper motion data. 
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Table 1 

Summary of the H30a Line Observations at v = 231.901 GHz 



Date 


t obs on Sgr A* 
(hr) 


T sy3 
(K) 


Array 3, 




Calibrators 


2006-04-08 


1.79 


96 


C 


3C273 b , 


Callisto c 


2006-05-28 


3.99 


173 


E 


3C279 b , 


Uranus b , J1924-292 b , Callisto b ' c 


2006-07-17 


2.19 


146 


V 


3C279 b , 


J1924-292 b , Callisto c 


2007-03-31 


3.67 


108 


C 


3C273 b , 


Callisto c 


2007-04-01 


3.35 


126 


C 


3C273 b , 


Callisto c 


2007-04-03 


1.99 


93 


C 


3C273 b , 


Callisto c 


2007-04-04 


1.93 


103 


C 


3C273 b , 


Callisto c 


2007-08-14 


3.63 


115 


E 


3C279 b , 


NRAO 530 c 


2007-08-20 


4.94 


95 


E 


3C454.3 


b , NRAO 530 c 


2007-08-31 


4.19 


134 


E 


3C454.3 


b , NRAO 530° 


2008-05-05 


2.67 


156 


C 


3C273 b , Ganymede" 



a Array configuration: C — Compact; E — Extended; V — Very Extended 

b Calibrator: bandpass 

c Calibrator: flux density scale 
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Table 2 

The H30a & H92a Line Measurements 



Radio Properties in IRS Regions 

_ . . a c c c Sh30c, Sh92c Sh30o a 

Region A« AO Sh30q S 2 2Gh z -3 -7; t; Vlsr AVfwhm 

_ »22GHa i'22GHz i'H92a 

(arcscc) (arcscc) (Jy beam ) (Jy beam ) (km s ) (km s ) 



IRS 1W 


5 


.27 


0.57 





21±0 


.008 


0. 


16±0. 


008 


1.3±0.08 


0.06±0.005 


21±1 


14±1 


46±2 


IRS 10W 


6 


..38 


5.02 





.15±0 


.018 


0. 


10±0. 


003 


1.5±0.19 


0.07±0.004 


21±1 


69±2 


39±5 


IRS 5 


8 


.53 


9.63 





10±0 


.019 


0. 


04±0. 


006 


2.5±0.6 


0.1±0.02 


25±2 


115±3 


33±8 


IRS 16 


2 


.30 


-0.26 





03±0 


.005 


0. 


08±0. 


015 


0.4±0.1 


0.026±0.008 


15±4 


-69±10 


141±27 


IRS 6 


-6 


.37 


1.00 





,12±0 


.011 


0. 


15±0. 


006 


0.80±0.08 


0.034±0.003 


23±2 


-146±5 


109±13 


IRS13E 


-3 


.17 


-1.53 





,08±0 


.006 


0. 


25±0. 


030 


0.32±0.05 


0.010±0.002 


32±4 


-31±3 


86±8 


IRS 2L 


-3 


.57 


-3.93 





21±0 


.011 


0. 


23±0. 


030 


0.91±0.13 


0.031±0.005 


29±1 


-269±1 


55±3 


IRS 12N 


-3 


.50 


-6.80 





,04±0 


.010 


0. 


04±0. 


002 


0.76±0.22 


0.036±0.013 


21±6 


-229±7 


40±17 


IRS 33 





.30 


-2.65 





,04±0 


.008 


0. 


19±0. 


024 


0.21±0.05 


0.009±0.003 


24±5 


-195±12 


126±29 


IRS 21 


2 


.18 


-2.83 





,08±0 


.009 


0. 


16±0. 


018 


0.50±0.08 


0.023±0.005 


22±3 


-99±4 


80±11 


IRS 9W 


2 


.85 


-5.62 





,06±0 


.007 


0. 


11±0. 


011 


0.57±0.09 


b 


b 


288±4 


60±8 


IRS 9 


5 


.20 


-7.30 





,08±0 


.009 


0. 


09±0. 


009 


0.85±0.13 


0.043±0.008 


20±2 


174±5 


84±11 


IRS 4 


10 


.33 


-5.98 





,07±0 


.012 


0. 


05±0. 


006 


1.4±0.3 


0.07±0.02 


20±3 


164±4 


47±10 


IRS 28 


10 


.60 


-6.80 





10±0 


.010 





07±0. 


006 


1.5±0.2 


0.06±0.01 


25±2 


168±4 


49±22 


IRS 7 





03 


5.52 





03±0 


.010 


0. 


15±0. 


001 


2.0±0.6 


0.07±0.03 


28±8 


-124±10 


29±22 


Minicavity s 


-1 


.75 


-2.80 





,08±0 


.006 


0. 


18±0. 


020 


0.45±0.07 


0.016±0.002 


29±2 


6±2 


41±4 



Radio Properties of Sgr A West 



Component J SmoadV / Sn92adV S22GHZ Area 

(Jy km s _1 ) (Jy km s _1 ) (Jy) (arcmin 2 ) (pc 2 ) 



Sgr A West tot al c 1700±300 d 60±7 16.3±0.2 1.0±0.1 5.5±0.5 

A c 300±20 19±2 7.1±0.1 0.10±0.01 0.55±0.05 

B f 1400±300 41±7 9.2±0.2 1.0±0.1 5.5±0.5 

a Dctcrmincd from the H30q: line. 

b No H92a data at the line velocity. 

c Thc overall Sgr A West region (80 // x45 // ). 

d The line flux determined from the IRAM 30m observation. 

e The bright H30a line emission region (the Northern and Eastern Arms in Fig. 2). 
f The residual of Sgr A West after taking out the contribution from component A. 

s The line and continuum intensities listed in this row reflect averaged values of the emission from both the rim and the void of the minicavity 
in the 2 beam that is similar to the diameter of the minicavity. 
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Table 3 

Physical Conditions of Ionized Gas in IRS Regions 



Region 


T 
(10 3 K) 


n c 

(10 4 cm- 3 ) 


EM 
(10 7 cm" 6 pc) 


t c (22GHz) 

(io- 3 ) 


m' 


( ^ ) 


/ A yH92 Q \ * 


V A Vfwhm / 


V *Vg* ) 


IRS 1W 


6.7±0.5 


6.0±0.4 


1.8±0.1 


14. ±1.1 


0.71 


0.88 


0.82 


IRS 10W 


6.8±0.6 


5.3±0.4 


1.1±0.1 


8.5±1.0 


0.71 


0.85 


0.72 


IRS 5 


5.0±1.1 


7.2±1.2 


0.4±0.1 


4.7±1.3 


0.68 


0.78 


1.2 


IRS 16 


6.7±2.0 


2.7±0.7 


0.8±0.2 


6.4±2.7 


0.68 


0.99 


0.37 


IRS 6 


5.0±0.6 


11. ±1.1 


1.6±0.1 


18. ±2. 3 


0.71 


0.96 


1.8 


IRS 13E 


12. ±2. 5 


21. ±3. 6 


3.2±0.5 


12. ±1.9 


0.84 


0.84 


2.0 


IRS 2L 


8.5±1.3 


14. ±1.4 


2.7±0.4 


16. ±2. 9 


0.78 


0.80 


1.7 


IRS 12N 


12. ±4. 7 


2.8±1.0 


0.5±0.2 


1.8±0.6 


0.77 


0.80 


0.27 


IRS 33 


13. ±3. 3 


10.±3.0 


2.6±0.3 


8.5±1.8 


0.81 


0.97 


0.92 


IRS 21 


10. ±1.8 


8.5±1.3 


2.0±0.3 


9.1±2.0 


0.78 


0.94 


0.92 


IRS 9W 


11.-12. 


0.8-18 


1.3-1.3 


5.5-4.9 


0.74-0.83 


0.76-0.93 


1.7-0.08 


IRS 9 


5.7±0.9 


6.4±0.8 


1.0±0.2 


9.7±1.8 


0.69 


0.97 


0.97 


IRS 4 


5.8±1.3 


5.0±1.0 


0.5±0.1 


4.8±1.4 


0.68 


0.91 


0.75 


IRS 28 


5.6±0.8 


8.9±1.0 


0.7±0.1 


6.8±1.1 


0.71 


0.86 


1.4 


IRS 7 


7.0±2.5 


8.8±2.7 


0.2±0.1 


1.2±0.5 


0.74 


0.64 


1.2 


Minicavity 


9.3±0.9 


8.5±0.7 


2.0±0.3 


10.±2.0 


0.77 


0.78 


0.96 



To T 1 10.87 

^ — — ~ I 630(1 — — /33oT"c(232 GHz) sa 6 30 for a very small value of rc(232 GHz) in a region at the Galactic center. 



T* L 2 



*AVrJ H = J 81n ( 2 ) fcTe is thc thermal Dopplcr line width. 
M H 
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